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We consider different patterns of supersymmetry breaking gaugino masses, and implications of 
these patterns for the phenomenology of neutralinos and charginos in models of low energy su- 
persymmetry. We outline a general procedure for obtaining approximate values for the neutralino 
masses relevant for our analysis, and describe the constraints on the gaugino mass parameters which 
follow from the present experimental limits on the mass of the lightest chargino. We evaluate an 
upper bound on the mass of the lightest neutralino that follows from the structure of the mass 
matrix in different models for the gaugino mass parameters. Using the experimental lower limit for 
the chargino mass, we examine the lower bound for the neutralino masses. Using a sum rule for the 
squared masses of charginos and neutralinos, we discuss how to distinguish between different models 
for the supersymmetry breaking scenarios. We discuss in detail the decay modes of neutralinos and 
charginos in different models of supersymmetry breaking. Our analysis shows that by measuring 
the masses and decay properties of the neutralinos, one can distinguish between different patterns 
of supersymmetry breaking in the gaugino sector. We then compare the dark matter characteristics 
that arise in different models of supersymmetry breaking. 

PACS numbers: 12.60.Jv, 14.80.Nb, 14.80.Ly 
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I. INTRODUCTION 

Supersymmetry is at present a leading candidate for physics beyond the standard model (SM). In supersymmetric 
models the Higgs sector of the standard model, which is crucial for its internal consistency, is technically natural, 
which makes these models very appealing. Since in nature there are no supersymmetric particles with the same mass 
as ordinary particles, supersymmetry must be a broken symmetry at low energies. The specific mechanism which 
breaks supersymmetry is important in determining the masses of superpartners of the SM particles, and, hence, 
the experimental signatures of supersymmetry. It is widely expected that at least some of the superpartners will be 
produced at the Large Hadron Collider (LHC), which has started its operations. However, most of the supersymmetric 
particles that are likely to be produced at the LHC will not be detected as such, since they will eventually decay 
into the lightest supersymmetric particle (LSP), which is stable as long as the R-parity (i?p) is conserved. Thus, the 
experimental study of supersymmetry involves the study of cascade decays of the supersymmetric particles to the LSP, 
and the subsequent reconstruction of the decay chains. The LSP in a large class of supersymmetry breaking models 
is the lightest neutralino, which has, therefore, been a subject of intense study for a long time [J-lHl- A stable lightest 
neutralino is also an excellent candidate for dark matter Q. In view of the possible production of supersymmetric 
particles and their subsequent decay into the lightest neutralino at the LHC, the properties of the lightest neutralino, 
and also those of heavier neutralinos and charginos, which often appear in the cascade decays, are of considerable 
importance. In particular a detailed study of the lightest neutralino, especially the predictions for its mass, are of 
great importance for the supersymmetric phenomenology. 



* Electronic address: katri.huituOhelsinki.fi 

^ Electronic address: j.laamanen@science.ru.nl 

Electronic address: ppanditaQnehu. ac . in 

^ Electronic address: paavo.tiitolaOhelsinki.fi 



2 



In the minimal supersymmetric extension of the standard model (MSSM) at least two Higgs doublets Hi and H2 
with hypercharge (Y) having values —1 and +1, respectively, are required to generate masses for all the SM fermions 
and gauge bosons, and to cancel triangle anomalies. By the minimal supersymmetric extension of the SM we here 
mean the model with minimal particle content and the standard model gauge group. The fermionic partners of these 
Higgs doublets mix with the fermionic partners of the gauge bosons to produce four neutralino states Xi,i = 1, 2, 3^4, 
and two chargino states xt^"^ = ^i^. In extended supersymmetric models, there can be extra neutralino states 

The masses of the neutralinos and charginos depend, besides other model parameters, on the soft supersymmetry 
breaking gaugino masses corresponding to the SU{2)l and U{1)y gauge groups. At present there are several models 
of supersymmetry breaking. The model of supersymmetry breaking that has been studied most extensively is the 
gravity mediated [loj supersymmetry breaking model. In this class of models, supersymmetry is assumed to be broken 
in a hidden sector by fields which interact with the SM particles and their superpartners (the visible sector) via only 
the gravitational interactions. Whereas this mechanism of supersymmetry breaking is simple and appealing, it suffers 
from the supersymmetric flavor problem. On the other hand, in a different class of models supersymmetry is 
broken in a hidden sector and transmitted to the visible sector via SM gauge interactions of messenger particles. This 
mechanism of supersymmetry breaking, the gauge mediated supersymmetry breaking, provides an appealing solution 
to the supersymmetric flavor problem. Both, the gravity mediated and gauge mediated supersymmetry breaking 
(GMSB) models, have their distinct experimental signatures. In the GMSB models, the lightest supersymmetric 
particle is gravitino instead of a neutralino. 

The soft supersymmetry breaking terms in the two breaking mechanisms described above have contributions origi- 
nating from the super- Weyl anomaly via loop effects. If gravity and gauge mediation of supersymmetry breakin g ar e 
somehow suppressed, the anomaly mediated contributions can dominate, as may happen, e.g., in brane models |12| . 




If this happens, then this mechanism of supersymmetry breaking is referred to as anomaly mediated supersymmetry 
breaking (AMSB). Anomaly mediation is a predictive framework for supersymmetry breaking in which the breaking 
of scale invariance mediates between hidden and visible sectors. 

There is another simple pattern of supersymmetry breaking, which is a hybrid between anomaly mediated supersym- 
metry breaking and mSUGRA pattern, that arises from the mirage (or mixed modulus) mediation of supersymmetry 
breaking, which has low energy values for masses quite distinct from either of the two. Mirage mediation is naturally 
realized in KKLT-type moduli stabilization [l3| and its generalizations, a well known example being KKLT moduli 
stabilization in type IIB string theory [l^ . 

The soft gaugino masses provide a handle for identifying the type of supersymmetry breaking [T5l - [l7| in the gaugino 
sector. In [l^ explicit examples of models of each of the patterns mentioned above were discussed. The possibility 
to detect the gaugino mass nonuniversality at the LHC was studied in [l^. It is because of the distinctive patterns 
of gaugino masses that one is tempted to believe that neutralinos and charginos are a key in understanding the 
supersymmetry breaking mechanism. 

In this paper we consider the phenomenology of neutralinos and charginos in different models for the soft supersym- 
metry breaking gaugino masses. This includes the gravity mediated supersymmetry breaking, the anomaly mediated 
supersymmetry breaking, and the mirage mediation of supersymmetry breaking. In Section II we recall the essential 
features of the neutralino mixing and the resulting mass matrix, and the constraints on the parameters of the mass 
matrix that follow from the experimental lower bound on the mass of the lightest chargino. We then discuss the pat- 
terns of the gaugino masses that arise in the different patterns of supersymmetry breaking in the gaugino sector, and 
compute the neutralino mass spectrum and discuss the distinguishing features of the spectrum in different patterns of 
supersymmetry breaking gaugino masses. In Section III we discuss a general upper bound on the mass of the lightest 
neutralino, and evaluate this bound for different models for the soft gaugino masses. Here we also discuss an upper 
bound on the lightest neutralino mass that follows from radiative electrowcak symmetry breaking, which determines 
the absolute value of the Higgs(ino) mass parameter /x. In Section IV we discuss and evaluate a sum rule involving the 
squared masses of the charginos and neutralinos, and describe how this sum rule can be used to distinguish between 
different models of supersymmetry breaking gaugino masses. In Section V we study the two-body decays of charginos 
and neutralinos in different models, and discuss how the branching ratios for these decays can be used to distinguish 
between the underlying pattern of supersymmetry breaking gaugino masses. In Section VI we compute and compare 
the relic abundance of the lightest neutralinos in different models assuming that the lightest neutralino is the lightest 
supersymmetric particle. We discuss the constraints imposed on the parameter space by the precise limits on the relic 
density obtained by the Wilkinson Microwave Anisotropy Probe (WMAP) satellite [1^. Finally, in Section VII we 
summarize our conclusions. 
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II. CHARGINO AND NEUTRALINO MASSES AND GAUGINO MASS PATTERNS 

111 this section we will describe constraints on the parameters of the neutralino mass matrix which follow from the 
present experimental limits on the mass of the lightest chargino. We will then discuss the patterns for the gaugino 
mass parameters that arise in different supersymmetry breaking scenarios, and the resulting consequences for the mass 
of the lightest neutralino. For completeness, we have summarized the neutralino and chargino mass matrices [l9j in 
the Appendix, and outlined the procedure to obtain approximate eigenvalues for the neutralino mass matrix. 

A. Experimental Constraints 

Collider experiments have searched for the supersymmetric partners of the standard model particles. No super- 
symmetric partners of the SM particles have been found in these experiments. At present only lower limits on their 
masses have have been obtained. In particular, the search for the lightest chargino state at LEP have yielded lower 
limits on its mass [1^. The lower limit depends on the spectrum of the model [2l[. Assuming that toq is large, the 
limit on the lightest chargino mass following from nonobservation of chargino pair production in e'^ e'~ collisions is 

Af-± > 103 GcV. (2.1) 

The bound depends on the sneutrino mass. For a sneutrino mass below 200 GeV, the bound becomes weaker, since 
the production of a chargino pair becomes more rare due to the negative interference between 7 or Z in the s-channel 
and v in the t-channel. In the models we consider, m,> is close to toq- When < 200 GeV, but to^ > m^±, the 
limit becomes [2lj 

A/-± > 85 GeV. (2.2) 

For the parameters of the chargino mass matrix the limit (j2.1[) implies an approximate lower limit [2^ [23j 

M2, > 100 GcV. (2.3) 

The limits Eq. (|2.3p on the parameters A'/2 and fi arc found from scanning over the MSSM parameter space and are 
thus model independent. 

B. Gaugino Mass Patterns 

Having constrained the parameters M2 and ^, which enter the chargino as well as the neutralino mass matrix, 
we now turn to the theoretical models for the supersymmetry breaking gaugino mass parameters Mi,M2, and M3. 
Theoretically, a simple set of patterns has emerged for these SUSY breaking parameters, which can be described as 
follows. 

1. Gravity mediated supersymmetry breaking 

The first pattern, which has been the object of extensive studies, is the one which arises in the gravity mediated 
supersymmetry breaking models We recall that the soft gaugino masses Mi and the gauge couplings gi satisfy 
the renormalization group equations (RGE's) {\M^\ = Mg, the gluino mass) 

odM, , /33 \ , , 

Wn^^ ^ 2b,AUgl b,= i-,l,-3j, (2.4) 

16^'^ = hgf, t = ln(M/Mo), (2.5) 

at the leading order, where i = 1, 2, 3 refer to the C/(l)y, SU{2)]^ and the SU{3) gauge groups, respectively, and /i is 

the renormalization scale with /iq as a reference scale. Furthermore, gi = \J^9\ 92 = 9, and 53 is the SU{3)c gauge 
coupling. In the minimal supersymmetric standard model with gravity mediated supersymmetry breaking and with 
a universal gaugino mass mi/2 at the grand unified scale (GUT), usually referred to as mSUGRA scenario, we have 
the boundary conditions (a^ = gf /Air, i = 1, 2, 3) 

Ml = Af2 = A/3 = mi/2, (2.6) 
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Oil = Q!2 = Q!3 = OLG, (2.7) 

at the GUT scale Mq. The RGB's (|2.4p and (|2.5|) imply that the soft supersymmetry breaking gaugino masses scale 
like gauge couplings: 

ai(Mz) a2(Afz) aaiMz) aa ' ^ ' ' 

which implies that Mi / gf docs not run at the one- loop level. Although in the context of the gravity mediated 
supersymmetry breaking models arbitrary soft gaugino masses are possible, we shall here consider the mSUGRA 
realisation (|2.6|) of the gravity mediated supersymmetry breaking scheme. The relation p.8|) reduces the three 
gaugino mass parameters to one, which we take to be the gluino mass Mg. The other gaugino mass parameters are 
then determined through 



5q; 

Mi{Mz) - Mg - 0.14 Afg, (2.9) 

Sas cos^ tlw 



AUMz) = r^— Mg ~ 0.28 Mg, (2.10) 

as sm Ow 

where we have used the value of various couplings at the Z° mass 

a-\Mz) = 127.9, sin^ % = 0.23, a^iMz) = 0.12. (2.11) 

For the gaugino mass parameters this leads to the ratio 

Ml : A/2 : A^s ^ 1 : 2 : 7.1. (2.12) 

This pattern is typical of any scheme obeying Eqs. (|2.4p and (|2.6p . Note that the masses above are the running masses 
evaluated at the electroweak scale, Mz- This discussion of the gaugino mass parameters is valid at tree level. When 
the radiative corrections are included, the ratio for these parameters in mSUGRA is modified to 

Ml : A/2 : A/3 ~ 1 : 1.9 : 6.2. (2.13) 

Using the ratio (|2.13p and the lower limit (|2.3p . we have the constraint 

Ml > 50 GeV, (2.14) 

in the gravity mediated supersymmetry breaking models. 

It is important to point out here that whereas the mechanism of gravity mediated supersymmetry breaking is 
simple and appealing, it suffers from the supersymmetric flavor problem. On the other hand in a different class of 
models, supersymmetry is broken in the hidden sector and transmitted to the visible sector via Standard Model gauge 
interactions of messanger particles. This mechanism of supersymmetry breaking, the gauge mediated supersymmetry 



breaking ll|, provides an appealing solution to supersymmetric flavor problem. In order to maintain the successful 
gauge coupling unification of the minimal supersymmetric standard model, it is usually assumed that the messanger 
particles form a full GUT multiplet, in which case the resulting gaugino masses follow the mSUGRA pattern (|2.12p 
as a consequence of the assumption of gauge coupling unification at the GUT scale. However, there are more general 
gauge mediated supersymmetry breaking models [2J] which allow any hierarchy of gaugino masses. However, we will 
confine here to the minimal models in which the gaugino masses follow the mSUGRA pattern (|2.12p . 



2. Anomaly mediated supersymmetry breaking 

The second pattern of gaugino masses, which is distinct from the mSUGRA pattern and emerges under theoretical 
assumptions that arc appealing, arises in the anomaly mediated supersymmetry breaking models. Since the soft 
supersymmetry breaking parameters are determined by the breaking of the scale invariance, they can be written in 
terms of the beta functions and anomalous dimensions in the form of relations which hold at all energies. In the 
minimal supersymmetric standard model (MSSM), the pure anomaly mediated contributions to the soft supersym- 
metry breaking parameters Mx (gaugino mass), mf (soft scalar mass squared), and Ay (the trilinear supersymmetry 
breaking coupling, where y refers to the Yukawa coupling) can be written as 



A/a = —m3/2, 
9 



(2.15) 
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mj = -\{^^+'-^Py)ml,,^ (2-16) 
Ay = -^^3/2, (2.17) 

where is the gravitino mass, /3's are the relevant j3 functions, and 7's are the anomalous dimensions of the 

corresponding chiral superfields. An immediate consequence of these relations is that supersymmetry breaking terms 
are completely insensitive to physics in the ultraviolet. The degrees of freedom that are excited at a given energy 
determine the anomalous dimensions and beta functions, thus completely specifying the soft supersymmetry breaking 
parameters at that energy. We note that the gaugino masses are proportional to their corresponding gauge group 
P functions with the lightest supersymmctric particle being mainly a wino. However, it turns out that the pure 
scalar mass-squared anomaly contribution for sleptons is negative [25| . There are a number of proposals for resolving 
this problem of tachyonic slepton masses (26l - l30| , but some of the solutions may spoil the most attractive feature of 
the anomaly mediated models, i.e., the renormalization group (RG) invariance of the soft terms and the consequent 
ultraviolet inscnsitivity of the mass spectrum. A simple phcnomcnologically attractive way of parametrizing the 
nonanomaly mediated contributions to the slepton masses, so as to cure their tachyonic spectrum, is to add a common 
mass parameter toq to all the squared scalar masses [3l|, assuming that such an addition does not reintroduce the 
supersymmctric flavor problem. Such an addition of a nonanomaly mediated term destroys the attractive feature of 
the RG invariance of soft masses. However, the RG evolution of the resulting model, nevertheless, inherits some of 
the simplicity of the pure anomaly mediated relations. 

There are several alternative ways to generate these extra contributions to the soft squared masses in the anomaly 
mediated supersymmetry breaking scheme. In particular there are models of supersymmetry breaking mediated 
through a small extra dimension, where SM matter multiplets and a supersymmetry breaking hidden sector are 
confined to opposite four-dimensional boundaries while gauge multiplets lie in the bulk. We note that in this scenario 
the soft gaugino mass terms are due to the anomaly mediated supersymmetry breaking, and, therefore, are governed 
by (|2.15p . On the other hand, scalar masses get contributions from both anomaly mediation and a tiny hard breaking 
of supersymmetry by operators on the hidden sector boundary. These operators contribute to scalar masses at one 
loop and this contribution is dominant, thereby making all squared scalar masses positive. The gaugino spectrum is 
unaltered, and the model resembles an anomaly mediated supersymmetry breaking model with nonuniversal scalar 
masses [33 |. 

Using Eq. (|2.15p . we then have the following pattern for the ratio of the gaugino masses at tree level: 

Ml ■■ M2 : M3 ^ 3.3 : 1 : 9, (2.18) 

which, after radiative corrections (assuming 7713/2 ~ 40 TeV) are included, becomes 

Ml : M2 : Ma ~ 2.8 : 1 : 7.1, (2.19) 

in the minimal supersymmctric standard model with anomaly mediated supersymmetry breaking. Schemes in which 
this pattern is realized require a strict separation of hidden sector that breaks SUSY from the visible sector of the 
MSSM. This implies a strong sequestering, and requires that all supersymmetry breaking fields are sequestered from 
the visible sector. Nevertheless, it may be achieved in certain class of theories with extra dimensions or a conformal 
field theory sector. 

Using (|2.3p . and the anomaly pattern of the gaugino masses (|2.19[) . we have 

Ml > 280 GeV. (2.20) 
This is to be contrasted with the corresponding result (j2.14p for the gravity mediated supersymmetry breaking. 



3. Mirage mediated supersymmetry breaking 

There is a third simple gaugino mass pattern that arises from the mirage (or mixed modulus) mediation super- 
symmetry breaking, which is a hybrid between anomaly mediated supersymmetry breaking and mSUGRA pattern, 
and has low energy values for masses quite distinct from either of the two. Mirage mediation is naturally realized 
in KKLT-type moduli stabilization and its generalizations, a well known example being KKLT moduli stabilization 
in type HB string theory [l3|. Phenomenology and cosmology of mirage mediation have been studied in l33l-4ll. 
Signatures of the scenario at LHC and the spectrum of neutralino mass in particular have been studied in |15l . |42| . 
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The boundary conditions for the soft supersymmetry breaking terms that produce the mirage mediation scheme can 
be written as 



Ma - Ah 



hl(Mp//m3/2) 2 

^ + 16^^ ^''^'^ 

[ai + aj + flfc) j^-^ (7, + 7j + 7^)0; 



ln(A/p;/m3/2) 



hl(A/p;/?7l3/2) \ . 2 

16^^ 



(2.21) 



where Mq ^ 1 TeV is a mass parameter characterizing the moduh mediation, Mpi is the reduced Planck mass, ga are 
the gauge couplings and the corresponding one-loop beta function coefficients, 7^ arc the anomalous dimensions, 
and a = m^/2/[Mol'n{Mpi/m^/2)] = C(l) is a parameter representing the ratio of anomaly mediation to moduli 
mediation, and 



a jk 

2 d-ji 



d log /i ' 



(2.22) 
(2.23) 



where C2{4>i) is the quadratic Casimir operator for the gauge group with gauge coupling ga-, and has a value {N'^ 
l)/(2Ar) for the SU{N) gauge group. 

Thus the generic mirage mediation is parametrized by 



Mo, a, a;, Ci = 1 - Ui, tan/3. 



(2.24) 



The parameter values = = 1 and a = 1 correspond to the minimal KKLT compactification of type IIB theory 
with modular weight = 0, but other parameter values are also possible for different scenarios, for example choice 
of a = 2 with an^ — <^Hu — ^-nd ajj^ + aq^ = cu^ + cq^ = can possibly alleviate the fine tuning problem for the 
electroweak symmetry breaking [iJl. In our studies we have used the values Ci — — \. At low energies, the gaugino 
masses in mirage mediation can be written as 



Mjjj)_ 



1 



\n{Mpi/m3/2) 2 , \ Mo 
dGurOaa ) -2 



167r2 



9gut 

This leads to a unification of the soft gaugino masses at the mirage messenger scale (4£ 

q/2 



(2.25) 



M„ 



M, 



GUT 



^3/2 
Mpi 



(2.26) 



which is lower than GUT scale for positive values of a. For g^i/rp ~ 1/2 the resulting low energy values yield the 
mirage mass pattern 



Ml : M2 : M3 ~ (1 + 0.66a) : (2 + 0.2a) : (6 - 1.8a). 

When the radiative corrections are included for the supersymmetry breaking gaugino masses, we obtain 

Ml : M2 : Af3 ~ 1 : 1.5 : 2.1 for a = 1, 
Ml : M2 : Af3 ~ 1 : 1.2 : 0.92 for a = 2. 

where we have used the value Mq = 1 TeV. Thus, for the mirage mediation, we find 

Ml > 67 GeV for a = 1, 
Ml > 83 GeV for a = 2. 



(2.27) 



(2.28) 
(2.29) 



(2.30) 
(2.31) 
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Parameters 


mSUGRA 


AMSB 


Mirage a = 1 


Mirage a = 2 


tan 13 = 5, mo = 200 GeV 
tan 13 = 20, mo = 200 GeV 
tan ^ = 5, mo = 1 TeV 
tan /3 = 20, mo = 1 TeV 


(58,105,250,277) 

(103,278) 
(58,104,229,253) 

(103,253) 
(55,103,346,363) 

(103, 365) 
(58,104,211,240) 

(103,242) 


(85,245,505,518) 

(85,518) 
(85,237,474,482) 

(85,484) 
(102,286,629,638) 

(103,640) 
(103,286,534,541) 

(103,545) 


(72,85,165,176) 
(85,184) 

(72,94,173,197) 
(85,202) 


(163,186,473,489) 
(174,479) 

(140,161,549,566) 
(150,553) 



TABLE I: The lower limits on the masses of the four neutralino states and two chargino states in each model in the form 
(m^o , m^o , m^o, m-,^o) [GeV] followed by (m^±,m^±) [GeV], with the given set of parameters, following from the experimental 
lower bound on the mass of the lightest chargino. For the mirage mediation model with a = 2 the limit is not from the chargino 
mass bound, but from the requirement of the nontachyonic spectrum. 



4-. Comparison of the patterns of neutralino and chargino masses 

In Table |l] we show the lightest neutralino and chargino masses, which satisfy the experimental limit for the 
mass of the lightest chargino 

m for a particular parameter point. These masses are calculated using SOFT- 
SUSY(v.3.0.13) [1^. The absolute value of the Higgsino mixing parameter is determined by the condition of radiative 
electroweak symmetry breaking (REWSB), and thus depends on the soft scalar mass parameter niQ. The parameter 
Too also enters the radiative corrections through the scalar masses. The masses of the neutralinos in TableUhave been 
calculated for toq ~ 200 GeV and niQ = 1 TeV for both the mSUGRA and AMSB models, and Ci — Ui = 1 for the 
mirage mediation. The sign of the fi parameter was chosen positive. Changing it to the negative value can reduce the 
neutralino masses by a few GeV's, but this may lead to conflict with the b —> s"f constraint. In mSUGRA the trilinear 
A-parameter was set to zero. Changing that to nonzero values also may decrease the lightest neutralino masses by a 
few GeV's. These masses demonstrate the effect of the sfermion spectrum on the neutralino and chargino masses for 
the models that we have studied in this paper. For anomaly mediated supersymmetry breaking, we see the familiar 
result that the lightest neutralino is closely degenerate with the lightest chargino. The neutralino spectrum in AMSB 
models is typically heavier than in the case of mSUGRA. In the case of mirage mediation, the mass difference of the 
lightest and heaviest neutralino masses is smaller as compared to this mass difference in the other models. In addition, 
especially in the a = 2 model, the /^-parameter is smaller as compared to its value in the other models, leading to 
larger mixing of the gaugino and Higgsino components. The much heavier spectrum in the a = 2 mirage mediation 
model is due to the tachyonic stops in the spectrum for the lighter particles. 

The masses of the neutralinos are plotted in Fig. [T]for the mSUGRA, AMSB and the mirage mediation scenarios, 
respectively. In mSUGRA the lightest neutralino is more than 80% bino, while the second lightest one is more than 
70% wino. In AMSB, Xi is almost 100% wino and X2 bino. In the mirage mediation pattern with a = 1, the 
compositions of the two lightest neutralinos are more evenly divided between bino and wino, but the lightest one 
is dominantly bino and the second lightest one wino. For small M2, also the Higgsino component is nonnegligible 
in both. For a — 2, the /i-parameter becomes relatively small, and both Xi and X2 are more than 90% Higgsinos. 
However, for the small values of M2, the LSP is not a neutralino as we will discuss later in Sec. VI (see Fig. 112b). 

The ratio of the mass parameters, IM3I/IM2I is known to differ drastically in different models. When the radiative 
corrections are taken into account, the ratios for each pattern discussed above is calculated to be 



IM3 



IM2 



3.3, 



mSUGRA 



I Mo 



= 7.1, 



AMSB 



I Mo 



1.4, 



IM3 



IM2 



0.77. 



(2.32) 



If one applies the ratio to the masses of the particles, one finds 



3.8-4.3, 



mSUGRA 



7.3-7.4, 



AMSB 



1.7, 



\mg\ 



0.9, 



(2.33) 



a=2 



with other parameter values as given in Table 1. Thus, very large mass ratio of gluino and chargino hints to an AMSB 
type breaking, while small value hints towards mirage type of supersymmetry breaking. 
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FIG. 1: Masses of the neutralinos in the (a) mSUGRA, (b) AMSB, and for mirage mediation models with (c) a = 1 and (d) 
Of = 2. Here tan/3 = 10, sgn(/i) = +1, Aa = 0, mo = 1 TeV (mSUGRA) and 1.5 TeV (AMSB). lep-denoted shading indicates 
the violation of the LEP sparticle mass limits. 



III. THE GENERAL UPPER BOUND ON THE MASS OF THE LIGHTEST NEUTRALINO 



In this Section we shall consider a general upper bound on the mass of the lightest neutralino that follows from the 
structure of the neutralino mass matrix. Since some of the neutralino masses resulting from diagonalization of the 
mass matrix (|A5P can be negative, we shall for our purposes consider the squared mass matrix M^M.. This squared 
mass matrix can be written as 



/ 



Ml 



Z '^^ 

Mid 



-MzSn,{MiCi3 + /iS^) MzCn,{M2Cp + ^1S p) 



-MzSw{MiCp + ^isp) 

MzCn,{M2Cp + ^1S|3) 



MzSw{Misp + fj.cp) \ 
-MzCw{M2Si3 + ^icp) 



(3.1) 



where cw = cos 9w , sw = sni9w,cp = cos/3 and sp = sin/3, respectively. An upper bound on the squared mass of 
the lightest neutralino Xi can be obtained by using the fact that the smallest eigenvalue of M^Mq is smaller than 
the smallest eigenvalue of its upper left 2x2 sub-matrix 



/ + M|sin2 



w 



Iz sm cos tiw 



(3.2) 



\ - Af| sin Qw cos Qw M| + Af| cos^ Bw 



thereby resulting in the tree-level upper bound [H 



M\ < - 
Xi - 2 



Ml 



2{Ml - Ml) Ml cos 2ew 



(3.3) 



9 




700 



m^± [GeV] 



(a) 



700 
600 
^ 500 
|< 400 
300 
200 
100 



100 



200 



300 



(b) 



400 



500 



600 



700 



m-± [GeV] 




300 400 500 
m^± [GeV] 

(c) 



700 




300 400 500 
m-± [GeV] 

(d) 



700 



FIG. 2: Upper limit (solid line) and the mass (dashed line) for the lightest neutralino mass for (a) mSUGRA with mo = 1 
TeV and A = 0, (b) AMSB with mo — 1 TeV, and mirage mediation with a — c — 1 and with (c) a = 1, and (d) with a — 2, 
calculated at tree-level (red/gray) and with radiative corrections added (black) as a function of the mass of the lighter chargino. 
In all the plots tan/3 — 10 and sign(/i) = +1. 



We emphasize that the upper bound (|3.3p is independent of the supersymmetry conserving parameter /i and also 
independent of tan/3, but depends on the supersymmetry breaking gaugino mass parameters Mi and M2. Despite 
this dependence on the unknown supersymmetry breaking parameters Eq. (|3.3p leads to a useful bound on M^o. 
An alternative bound on the mass of the lightest neutralino can be obtained by considering the bottom-right 2x2 
sub-matrix 



leading to an upper bound 



A/| cos^ + //^ -Af|cos,3sin/^ 
- A/| cos /3 sin /3 A/| sin^ l3 + fi"^ 



M^o < I/if 



This bound, when supplemented by the electroweak symmetry breaking condition 
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FIG. 3: The sum rule (|4.ip plotted for different gaugino mass patterns. 



leads to the upper bound 



{ml — 1712 tan^ (3) 1 



(tan^ 13-1) 



-Mi 



(3.7) 



The importance of this bound stems from the fact that it relates the two sectors, namely the supersymmetry breaking 
gaugino masses and the Higgs masses. One can then use the RG evolution for the parameters on the RHS of the 
bounds to evaluate the bound. Note that the RG equations for the parameters on the RHS of the bounds involve 
the gaugino masses, and will, therefore, depend on the boundary conditions for the gaugino masses, and hence on the 
different supersymmetry breaking models for the gaugino mass parameters. In Fig. [2]we have plotted the upper limits 
for the lightest neutralino mass following from ()3.3p for the different supersymmetry breaking models as a function 
of m ~± . For all the four models studied we plot the tree-level masses, and the two-loop radiatively corrected masses 

calculated using SOFTSUSY(v.3.0.13) Hi]. In Fig. [H, we plot also the upper hmit from Eq. dSJ]), which gives 
the lowest mass upper bound in this case. Note that the mass of the lighter of the charginos is close to the value 
of the /^-parameter. The neutralino masses have been calculated assuming tan/3 = 10, and the other parameters as 
indicated in the Figure. We note that all gaugino masses can receive radiative corrections up to 20%, and, thus, 
difference between tree level and radiatively corrected neutralino and chargino masses can be significant in all models 
studied in this paper, although the difference is not explicit in Fig. [5] due to similar magnitude of correction for both 
the neutralino and the chargino. 

IV. SUM RULES 

We recall that in the minimal AMSB model, the mass difference between the lightest chargino and the lightest 
neutralino is small. The close proximity of the lightest neutralino and chargino masses is a direct consequence of 
Eq. (|2.15p . which gives for the ratios of the gaugino mass parameters \Mi\ : \M2\ : \M^\ ~ 2.8 : 1 : 7.1 , after taking 



into account the next to leading order radiative corrections and the weak scale threshold corrections [31[ as in (|2.19p . 
Thus, the winos are the lightest neutralinos and charginos, and one would expect that the lightest chargino is only 
slightly heavier than the lightest neutralino in all anomaly mediated models. It is not feasible to obtain sum rules for 
the masses of the neutralino states, since the physical neutralino mass matrix is a 4 x 4 matrix. However, from the 
trace of the squares of the neutralino and chargino mass matrices, one obtains a sum rule, which does not contain the 
Higgs mixing parameter /i, but which is present in the mass matrices. The sum rule can be written as 



2 Vm 



2 



[Ml - Mf] + AM^ - 2M| 



(4.1) 



By using the gaugino mass pattern for a specific model the sum rule (|4.ip can be expressed as as function of any 
of the gaugino masses. This is shown in Fig. [3] The average mass difference in the AMSB is first positive, but then 
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quickly turns negative, while in the minimal SUGRA model it is always positive. In the mirage mediation model the 
behavior is determined by the parameter a with a low value leading to a mSUGRA-like curve. Increasing a decreases 
the gradient of the curve until a = 2.17 (the value for which Mi = AI2) leads to a constant positive value. We note 
that this sum rule could be used as a signature for different supersymmetry breaking models, and in the case of mirage 
mediation it might be useful for determining the value of a, which can be calculated from the sum rule for specific 
values of gaugino masses. 

V. DECAYS OF NEUTRALINOS AND CHARGINOS 

In this Section we discuss the decays of charginos and neutralinos in different supersymmetry breaking models 
that we have discussed in this paper. As noted earlier, charginos and neutralinos are mass eigenstates, which are 
mo del- dependent linear combinations of charged or neutral gauginos and Higgsinos. Since the mass matrices of 
charginos and neutralinos depend on parameters Mi and M2, which are model dependent, the decays will depend on 
the model under consideration. As such the decay patterns of charginos and neutralinos can serve as a window on 
the underlying supersymmetry breaking mechanism in the gaugino sector. Here we shall mostly focus on two-body 
tree-level decays of neutralinos and charginos, if they are kinematically possible. If the neutralino or the chargino 
is sufficiently heavy, then two-body decays into a IV, Z", or a Higgs boson and a lighter neutralino or chargino are 
the dominant decay modes. Since in supersymmetric models with minimal particle content, the lightest Higgs boson 
is relatively light, the two-body decay containing the light Higgs boson is expected to be the dominant decay mode 
over a large region of parameter space. However, if some squarks or sleptons are relatively light, the two body tree- 
level decays of a heavy neutralino or chargino to quark-squark or lepton-slepton may be important. However, these 
decays are phenomenologically less important at a hadron collider like LHC, where neutralinos and charginos would 
be produced from the decays of strongly interacting squarks and gluinos. Neutralinos and charginos, which are heavier 
than squarks, would be hard to study at a hadron collider. 

We recah that if a heavier x° (* = 2, 3, 4) or a chargino {j = 1, 2) is produced at a collider, it will decay via a 
cascade until the lightest neutralino ( x") is produced. Thus, we are here interested in the branching ratios for the 
two-body decays 

t^x° + z\ t^xf + w^, x+^x° + w^+, xt^x^ + z\ (5.1) 

Xl-^X°+Hl x"^xf+H^, xt^X° + H+, xt^xt+Hl (5.2) 

These two body decays will dominate any tree-level three-body decays mediated by virtual squark or slepton exchange. 
These decays will also dominate any two-body decay, which is forbidden at the tree level, but which can proceed via 
loops, such as — > Xj + 7- 

If some of the neutralinos and charginos are heavier than some of the squarks and sleptons, then the two-body 
decays 

X-^qq,ll (5.3) 

can compete with the two-body decays into W,Z'^,H discussed above. The analytical expressions for the branching 
ratios of charginos and neutralinos into W,Z'^, and Higgs bosons, as well as into squarks/sleptons for arbitrary 
neutralino and chargino mixing angles are given in (47| . 

For the evaluation of branching ratios of charginos and neutralinos, we have calculated the spectra of the super- 
symmetric particles using SOFTSUSY(v.3.0.13) l46ll . and the decays of the supersymmetric particles using SUSY- 
HIT(v.l.3 with SDECAY vl.Bb/HDECAY v3.4) [ii]. In the following analysis we have used the parameter values 
tan/3 = 10 and sign(/^) = +1 for all models. In addition, for the mSUGRA and AMSB models, we have used mo = 1 
TeV. For the mirage mediation of supersymmetry breaking, we have used a = 1/2 and a = 1 with Ci — ai ~ 1. The 
mirage mediation scenarios with higher value for a can lead to squark LSP or tachyons, as will be discussed in next 
Section. While calculating the decay rates of charginos and neutralinos, we have imposed the experimental constraints 
following from LEP sparticle mass limits and the LEP lower bound on the lightest Higgs boson mass. Thus, in each 
of the following figures, the green (dark in grayscale) color denotes the area where LEP sparticle mass limits are 
violated, and light blue color (light) denotes the area where the lightest Higgs mass is below 114 GeV. In addition, 
the vertical Isp-denoted line in the mirage a = 1 figures indicates the minimum mass for which the neutralino is the 
LSP. The bsg-denoted line represents the BR{b — > S7) constraint. In the anomaly mediated supersymmetry breaking 
models, the bsg line practically coincides with the Higgs mass limit for these parameters, and is not drawn explicitly 
in order to reduce too many curves in the figure. (In the mSUGRA figures the constraint is obeyed throughout the 
mass ranges.) The LEP and 6 — >■ S7 constraints are calculated with micrOmegas (see Sec. IVip . An asterisk after 
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FIG. 4: Branching ratios for the two body decays of ul to gauginos and Higgsinos in (a) mSUGRA model, (b) mirage mediation 
scenario with a — 0.5 and (c) a = 1, and (d) the AMSB scenario. Shadings specify the region where the LEP mass hmits are 
not satisfied (dark) or the lightest Higgs mass is below 114 GeV (light). See the text for more details. 



the decay mode indicates that the charge conjugated mode is also included in the plotted value. A label of the form 
Xixx/yy indicates that the plotted value includes both channels Xi^x and XiVUi ^^'^ that each channel contributes 
equally. 

We begin by commenting the production of charginos and neutralinos at the LHC. The direct pair production of 
these particles is rare [49| compared to the pair production of strongly interacting particles [soj . However, neutralinos 
and charginos are also produced in the cascade decays of squarks and giuinos. As an example, we study the decays 
of squarks to neutralinos and charginos. 

The decay branching ratio of the left-handed u-squark is depicted in Fig. |31 It is seen that if gluino remains 
lighter than the squark, the decay is dominantly via the strong coupling to gluino and quark, and then further gluino 
would decay to the lightest neutralino with a quark pair or a gluon. If the two-body decay to quark and gluino is 
kincniatically not possible, the squark decays mainly to the lighter chargino and quark, xt^, or to the second lightest 
neutralino and quark, X2^, in mSUGRA and mirage patterns. In the AMSB, the branching ratio to Xi"*^ is larger 
than to X2U This is due to the fact that in the mSUGRA and mirage patterns the neutralino X2 is mainly wino, 
while in the AlVISB pattern Xi is dominantly wino. The right-handed squark un decays most of the time to x^'" in 
the mSUGRA and mirage patterns, and to X2^ in the AlVISB case, since for uj^ the decay to the bino dominated 
neutralino is favoured. 

For the chosen set of parameters with tuq = 1 TeV, the squarks are relatively heavy, especially when the experi- 
mentally allowed range for 6 — > S7 is taken into account. Thus one might want to consider the third generation stops, 
which are lighter than the other squarks due to mixing. In Fig. [5] we have plotted the two-body decay modes of ti to 
neutralinos, charginos and giuinos as a function of the stop mass. One finds that especially in the case of mSUGRA 
and AMSB, the cascades produce heavier neutralinos and charginos, e.g. in the AMSB pattern more than 80 % of 
the lightest stop decays to or xfb. This is different from the two mirage mediation scenarios, where the lightness 
of stop causes it to decay directly to top and the neutralino LSP, or to b-quark, W-boson and neutralino through a 
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FIG. 5: Branching ratios for the largest decay modes of ii to gauginos and Higgsinos in (a) mSUGRA model, (b) mirage 
mediation scenario with a — 0.5 and (c) a = 1, and (d) the AMSB scenario. Shadings specify the region where the LEP mass 
limits are not satisfied (dark) or the lightest Higgs mass is below 114 GeV (light). See the text for more details. 



3-body process. 

In Fig. [6] we have plotted the dominant decay modes of X2 ^ function of its mass for mSUGRA, the mirage 
mediation models, and the AMSB, respectively. We note from Fig. [S] that the dominant mode of X2 mSUGRA is 
into the Xi a-nd the lightest Higgs boson, when kinematically possible, and into Xi ^^^d the Z° boson for a light X2- 
The same applies for the two mirage mediation scenarios; one can see the opening of the Higgs mode to be pushed 
gradually to higher neutralino masses as the mirage scenario further deviates from the pure mSUGRA by increasing 
a. With an increasing a, stop becomes lighter and can be the lightest supersymmetric particle. In the heavier end 
of the spectrum the lightness of stops enables the stop decay mode. In the AMSB models, where the light chargino 
and neutralino are almost degenerate in mass, the decay mode to the lighter chargino and W-boson dominates the 
decay to Xi a-^d the lightest Higgs boson, until the sfcrmion decay modes become kinematically possible and begin 
to dominate. 

In Fig. [7] we have plotted the dominant two body decay modes of the X3 a function of its mass for mSUGRA, 
the mirage mediation models, and the AMSB, respectively. For X3 the chargino-VF mode is available in mSUGRA, 
and it dominates. However, the X2^'^ mode starts competing with increasing mass of Xs- In the mirage mediation 
models the stop-top channel is available in addition to the XiW and the decay modes. Furthermore, in the mirage 
mediation scenario, the light stop-top channel takes the role of the dominant decay mode for a large mass range, 
where the two mirage scenarios give quite similar branching ratios for the X3 decay. This is in contrast to the AMSB 
scenario, where the XaW and modes dominate. 

In Fig. [8] we have plotted the dominant two body decay modes of the X4 as a function of mass of Xa ■ mSUGRA 
and AMSB scenarios the XiW mode dominates, whereas in the mirage mediation the stop-top channel has the largest 
branching fraction. For the relevant mass range, in the two mirage scenarios, the xi decays look similar. Furthermore, 
the light Higgs decay modes are present in each breaking scenario. 

In Fig. [HI the dominant decay modes of Xi^ are plotted as a function of its mass for the different supcrsymmetry 



14 



(a) 



3i 



0.6 



0.4 



0.2 



rrih < 114 GoV 
lep I 



Xiuu/ cc 



x°bb — ■ 



rir (*) 
Xidd/ss 



400 



800 

,0 [GcV] 



1200 



1600 




400 



800 

m~o [GcV] 



1200 




400 



800 1200 

rn-o [GcV] 



1600 



(d) 



t 



0.6 



0.4 



0.2 



mh<114GoV 

x> ■ 
xtWi.: 



e/?.e(*) 
St fx 
T2r (*) 



rir(*) -- -V - 
{iRfi(*) 



400 600 



800 1000 1200 

m^p [GcV] 



1400 1600 



FIG. 6: Branching ratios for the two body decays of X2 in (a-) mSUGRA model, (b) mirage mediation scenario with a — 0.5 
and (c) Q = 1, and (d) the AMSB scenario. Also the three body channels are shown where no two body decays are possible. 
The (*) in the channel label signifies that the channel includes also the charge conjugated final state. Shadings specify the 
region where the LEP mass limits are not satisfied (dark) or the lightest Higgs mass is below 114 GeV (light). 



breaking models. In mSUGRA the only possible two-body decay mode, until the slepton (stau) channel becomes 
available, is the xV^ mode. In the mirage mediation scenario the lightness of the stop enables the tih mode, which 
becomes the dominant one for a = 1. The other available channel is the XaW mode, which in fact is the only two-body 
mode in the a = 0.5 case, thus making the two mirage scenarios distinguishable with respect to the lighter chargino 
decay. Also, the neutralino LSP and the 6 — > 57 requirements push the chargino mass to be reasonably heavy in the 
a = 1 mirage scenario. In the AMSB models the lighter chargino mass is so close to the lightest neutralino mass that 
only three body decays are available. 

Finally, in Fig. [10] we plot the dominant two body decay modes of the X2 a function of its mass for the models 
discussed in this paper. In mSUGRA the W and Z'^ modes dominate for the low mass range, but eventually the 
light Higgs channel becomes a competing mode. In the AMSB scenario the light Higgs, and W modes have 
approximately the same branching ratio. In contrast to the other models, in the mirage mediation the tib mode 
dominates by a large margin. 

In general, it is obvious that in the AMSB models the close proximity of the lightest neutralino and the lighter 
chargino is important in the decays of heavier neutralinos and chargino as well. In addition to forcing the lighter 
chargino to decay via a three body mode, this feature makes the second lightest neutralino to decay dominantly into the 
lighter chargino and W boson or other charged sparticles. The predictions of the anomaly mediated supersymmetry 
breaking and the mSUGRA are very different here. 

In mirage mediation the mass spectrum for neutralinos and charginos is more tightly packed, and thus for X2 ^-nd xt^ 
it is expected that they decay via a three body mode, if the mass scale is light enough, which again is a distinguishing 
feature as compared to the mSUGRA and AMSB models. The lightness of stop allows in many cases the stop decay 
channel in the mirage mediation models, which then becomes the dominant decay mode. 
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FIG. 7: Branching ratios for the decays of xs in (a) mSUGRA model, (b) mirage mediation scenario with a — 0.5 (c) in the 
mirage mediation scenario with a = 1, and (d) the AMSB scenario. 

VI. RELIC DENSITY 

In this Section, we study the imphcations of different patterns of supersymmetry breaking gaugino masses for the 
reUc density of lightest neutralino, and the constraints imposed on the parameter space by the precise limits on 
the relic density obtained by the Wilkinson Microwave Anisotropy Probe (WMAP) satellite. Requiring the lightest 
neutralino to form all of the dark matter as a thermal relic is very limiting constraint on the parameter space, and 
it should be kept in mind that the possible dark matter might also be created non-thermally or the excess thermal 
production diluted for example by an entropy increase after the freeze-out [5ll - [53| . Therefore we refer to the WMAP 
constrained parameter space as a WMAP-preferred relic density area. 

The relic density in the mSUGRA [54| - |6G | and AMSB [3l|, l6ll - [64| scenarios has been studied extensively. Neutralino 
dark matter in the mirage mediation SUSY breaking model has been considered in j4ll . [65l - l67j . Here we consider the 
combined information from the relic density and decay modes of the second lightest neutralino. In Fig. [TT] the main 
decay modes of X2 ''^^^ the WMAP-preferred relic density stripe in the (M~±,tan/3) plane are plotted for mSUGRA 
(^0 = 0,mo = 120 GeV and mo = 1 TcV), AMSB (mo = 5 TeV), and mirage mediation (a = 1 with Ui = Ci = 1) 
for sign(/i) = -|-1. The spectrum is calculated using SOFTSUSY (v. 3. 0.13) [4Q|, and the relic density and constraints 
using micrOmegas (v.2.2.CPC.i) [g^. For the relic density, we use the WMAP combined three year limits [l^ 

f^cm//^' = 0.11054t;j°g^^^ (2a). (6.1) 

In the figures below, the filling denoted by wmap is the WMAP preferred region. For the b ^ sj experimental 
branching fraction, the two sigma world average has been used BR{b 57) = (355 ± 24t% ± 3) x 10"^, and for 
the particle masses the limits of [6^ are applied. In the figures, lep shows an area where the experimental sparticle 
mass limits are not met, rge shows an area where there are tachyons or no radiative EWSB, and Isp the area where 
neutralino is not the LSP. The curve mji = 114 GeV is depicted in the figures for the reference (dash-dotted line 
denoted by h). 
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FIG. 8: Branching ratios for the decays of xl in (a-) mSUGRA model, (b) mirage mediation scenario with a — 0.5 and (c) 
a = 1, and (d) the AMSB scenario. 



In the mSUGRA scenario the WMAP preferred relic density area can be obtained either for a smah mo with 
moderate tan /3, where the scalars are hght and stau coannihilates with hght neutrahno, or for a larger mo with quite 
large tan p, where Xi annihilates through Higgs channel. The lightest neutralino is very bino-like, which in general 
results in a large relic density. In Fig. [TTk the WMAP preferred region is shown for uiq = 120 GeV. This is achieved 
via eoannihilation with stau for 'm~± < 500 GeV and tan/3 < 20. The dominant decay mode along the WMAP 

stripe is to fir -pair until various other leptonie modes and the light Higgs mode take over with increasing xt mass. 
While the stau eoannihilation region exists also for the heavier bino masses, the WMAP preferred region becomes 
very narrow due to the relic density mass dependence on the LSP mass. Increasing mg to 200 GeV would raise the 
WMAP preferred region to around m-;± ~ 500 — 700 GeV and 12 < tan/3 < 25, but the stripe would be clearly 

Xi 

narrower due to stau eoannihilation with heavier bino. For mo = 1 TeV, Aq = 0, and /i > 0, the dark matter area 
can be found around tan/3 ~ 50 (t^], as shown in Fig. [TTb . Here both neutralino annihilation through a heavy Higgs 
resonance and stau eoannihilation are effective. Along the WMAP stripe Xi^i is the dominant decay mode until tit 
becomes kinematically accessible. 

We note that since in the AMSB scenario the lightest neutralino is almost a pure wino for a large region of the 
parameter space, the relic density for moderate mo values tends to be smaller than the WMAP observation. Therefore 
a spectrum with rather large LSP mass is required for the AMSB to provide the observed relic density. In Fig. [TTb 
the value of mo = 5 TeV has been used. The WMAP region occurs for m ± ~ 2 TeV, where the decay modes to 

Xi 

slepton-lepton pair dominate the X2"'^6cays. 

In mirage mediation scenario (Fig. Illd ) with a = 1 with the mentioned choice of parameters the lightest neutralino 
is mostly bino-like. Therefore the existence of the WMAP-preferred relic density region in general relies on the 
eoannihilation with the NSLP, which for the large part of the parameter space is stop. (For high tan P the NLSP 
near the Isp-denoted region is stau.) Taking into account the b ^ sj constraint and requiring that the relic density 
agrees with the WMAP observations practically removes the WMAP-preferred region. If a ten percent error in the 
theoretical calculation of the 6 — >■ 57 constraint is assumed, the WMAP-stripe on the right of the bsg th. -denoted line 
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FIG. 9: Branching ratios for the decays of xt in (s-) mSUGRA model, (b) mirage mediation scenario with a — 0.5 and (c) 
a = 1, and (d) the AMSB scenario. 



becomes allowed. The peaking behavior in the relic density is due to the efficient s-channel annihilation through heavy 
Higgs bosons. Around tan /3 ~ 10 the annihilation rate is enhanced via stops decaying through the Higgs resonance, 
while around tan/3 = 18 it is the lightest neutralino which annihilates directly though the resonance. For a large 
region of the parameter space the NLSP is the lighter stop, so the WlVIAP-preferred relic density parameter space 
follows the neutralino co-annihilation region with stop. The figure is divided into areas where the dominant X2 decay 
mode exceeds 50 %. If the constraints are obeyed, the dominant decay mode for X2 to Xil^ or possibly to tii{*) for 
the WlMAP-preferred relic density area. 

In Fig. [TTn . the composition of the lightest neutralino and the WlVIAP-prcfcrred relic density stripe in the mirage 
mediation scenario in the (71/3/2,0;) plane are plotted (for tan/3 — 10 and sign(/i) = +1) with = c; — 1. The 
lightest neutralino composition is painted in RGB-color encoding (i.e., colors (or hues of gray) indicate the particle 
as shown in the figure color bar; therefore the mixture of the colors (or the shading in black and white) describes the 
nature of the Xi-composition). As the a-parameter increases, the lightest neutralino composition makes a transition 
from the bino-like to very Higgsino-like. Above the Isp-denoted line the stop is the LSP, and therefore that area is 
not suitable for the dark matter considerations. However, at sufficiently high alpha value the REWSB forces the 
/i-parameter to very light value, hence reducing the Higgsino like neutralino mass below the sparticle masses and the 
lightest neutralino becomes the LSP again. The areas can be seen in Fig. [T^ . where also the domains of X2 decay 
modes exceeding 50 % are marked. The thin WMAP-preferred relic density band lies along the line where the lightest 
neutralino can coannihilate with the lighter stop, hence reducing the relic density to an acceptable level. Here the X2 
decays dominantly to Xi^ or to fir pair, as indicated in the figure. 
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FIG. 10: Branching ratios for the decays of xt in mSUGRA model, (b) mirage mediation scenario with a — 0.5 and (c) 
a = 1, and (d) the AMSB scenario. 



VII. SUMMARY AND DISCUSSION 



The lightest neutralino is expected to be the lightest supersymmetric particle in supersymmetric models with im- 
parity conservation. It is expected to be the end product of decays of supersymmetric partners of the Standard 
IVIodel particles that are likely to be produced at the Large Hadron Collider. Thus, its mass, and its properties are 
of considerable importance for the supersymmetric phenomenology. The mass of the lightest neutralino, as well as 
those of its heavier partners, depend on the mechanism of supersymmetry breaking in the gaugino sector. 

We have carried out a detailed study of the spectrum of neutralinos and charginos in different models of super- 
symmetry breaking gaugino masses. Because of its importance, we have investigated the properties of the lightest 
neutralino in different patterns of supersymmetry breaking in the gaugino sector. We have calculated lower limits on 
the masses of the neutralinos and the charginos, taking into account the current experimental limits on the mass of 
the lightest chargino. Although these limits depend, through radiative corrections, on parameters other than gaugino 
sector parameters, we have found that this dependence is mild, and thus the limits for the neutralino and chargino 
masses can be considered to be relatively robust. 

We have calculated an upper bound on the mass of the lightest neutralino as a function of the lightest chargino 
mass. We see that for the models of supersymmetry breaking considered in this paper, only in the mirage mediation 
model with large a, the upper bound found from the lower right hand two-by-two part of the mass matrix becomes 
relevant. 

The sum rule for the neutralino and chargino squared masses is one of the distinguishing features of the supersym- 
metry breaking mechanism in the gaugino sector. For mSUGRA and AlVISB the sign of the average squared mass 
difference for charginos and neutralinos is different for Afi larger than 0(100 GeV), and for mirage mediation the 
parameter a can be deduced from the sum rule. 

We have also discussed in detail the decay patterns of the neutralinos and charginos in different models. An interest- 
ing result of our work is that detection of neutralino and chargino decay patterns gives important information on the 



19 




1600 



1800 



2000 



2200 



2400 



2600 



400 



600 



800 



1000 



1200 1400 



m ± [GcV] 



m ± [GeV] 



FIG. 11: Contours of constant branching ratio for the leading two-body decay modes of X2 superposed on the same plot with 
several constraints for (a) mSUGRA scenario for mo = 120 GeV, (b) mo = 1 TeV, (c) AMSB for mo = 5 TeV and (d) the 
mirage mediation for a = 1. The — > 57 constraint is obeyed right of the dotted bsg-denoted line and the lightest Higgs mass 
is more than 114 GeV on the right of the h denoted dash-dotted line. 



nature of the underlying supersymmetry breaking mechanism, and may help in identifying the correct supersymmetry 
breaking pattern. 

In Section V it was shown that the second lightest neutralino and the lighter chargino arc produced in large amounts 
in squark decays. This is interesting, since a promising signal to detect weakly interacting particles at Tcvatron and 
at LHC is considered to be the associated production xfx2i ^^e e.g. [tiIIt^ and references therein. Let us consider 
produced xf, X2 view of the cascade decays in Figs. d!])-®. It is seen that in the studied breaking patterns the 
largeness of the trilepton signal varies significantly. In the mSUGRA pattern, ti decays to all the heavier neutralinos 
and charginos with nonnegligible branching fractions. The contribution ii — >■ Xi^^/X2^ is at a few percent level, but 
more events come from the decays of X3 4: xt ■ Thus from iiii production there is an additional contribution to the 
trilepton signal, accompanied by a number of jets. In the AMSB pattern, the enhancement of trileptons is significant, 
ti's decay 60% of the time to and 20 % of the time to xt^- soon as kinematically possible, the X2 decays to 
a slepton and lepton, and Xi decays leptonically 25% of the time. In mirage pattern, stops tend to decay directly to 
the lightest neutralino and no enhancement is expected. 

Since the lightest neutralino is a possible candidate for the particle dark matter, we have calculated its relic density 
in different supersymmetry breaking models combining the information coming from decay patterns. While in the 
mSUGRA model typically a narrow range with the observed relic density occurs, in the AMSB model the relic density 
remains below the WMAP limit for the sub- TeV scale spectrum. In mirage mediation models the observed dark 
matter range is narrow and close to the stop LSP region, unless the heavy Higgs resonance can be found. We note 
that it is not necessary that neutralino is the only dark matter particle, even if it were the lightest supersymmetric 
particle. Furthermore, it is possible that the R-parity is broken at least slightly in nature. This would lead to the 
neutralino decay, even if the breaking were so tiny that it would not show up in the experiments. 
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FIG. 12: Lightest neutralino composition (a) and the leading X2 decay modes (b) in the mirage mediation scenario in (M3/2, o) 
plane for sgn(/i) = +1, tan/3 = 10 and ai = Ci = 1. The narrow light yellow band in (a) (red in (b)) indicates the WMAP 
preferred relic density area. The b ^ s'y constraint is obeyed below the dotted bsg-denoted Une and the lightest Higgs mass 
is more than 114 GeV below the h denoted dash-dotted hne. The Isp denoted (light blue) line near the WMAP filling limits 
the area, above which the lightest neutralino is not the LSP except for the area near a = 2, which can better be seen in (b). 
The black areas limited by the rge-denoted line depicts the area where there are either tachyons or no REWSB. In (b) the 
domains of branching ratio exceeding 50 % for the leading decay modes of X2 B.ie drawn for the same parameters, including 
the constraints. 
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Appendix A: Chargino and Neutralino Mass Matrices 



In the wino-Higgsino basis 



(-iA+,VHj, V-/ = (-«A-,V^^J, J = 1,2, 



(Al) 



where = (1/V2)(\^ =F A^), and the superscripts 1, 2 refer to SU{2)l indices, the chargino mass matrix can be 
written as llOll 



M+ = 



Mo 



^/2Mws\x\l3 



V2Mw cos l3 



(A2) 



where M2 is the supersymmetry breaking SU{2)l gaugino mass, fi is the Higgs(ino) mixing parameter, and tan/3 is 
the ratio of the vacuum expectation values of the neutral components of the two Higgs doublets H2 and Hi . We shall 
denote the eigenstates of the chargino mass matrix (jA2p as xf ^.nd X2 ^ with eigenvalues M^± , respectively. The 

eigenvalues are most easily obtained from the diagonalization of M.±^ M.± resulting in the squares of the chargino 
masses 



Mi 



2mly f - 4{M2fi - sin 2^^ 



(A3) 
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On the other, in the bino-wino-Higgsino basis 

^- = (-«A', - ^A^ J = 1, 2, 3, 4, 



(A4) 



where A' and A'^ are the two-component gaugino states corresponding to the U{1)y and the third component of the 
SU{2)l gauge groups, respectively, and iphi^ ^If2 two-component Higgsino states, the neutralino mass matrix 

can be written as 1191 



( Ml 

Ah 

— Mz cos I3suy9w Mz cos 13 cos 9 w 

\ Mz sin /3 sin 6*^/ —Mz sin f3 cos 9w 



— Mz cos (3 sin 9 w Mz sin /3 sin 6'vi/ \ 

Mz cos/?cos6'vK — Mz sin /3 cos 

-fi 

-fi / 



(A5) 



All is the supersymmetry breaking U{1)y gaugino mass, and g' and g are the gauge couphngs associated with the 
U{1)y and the SU{2)l gauge groups, respectively, with tan6'n/ = g' /g, and A/| = +u|)/2. The neutralino 

mass matrix can be diagonalized by a unitary transformation N 



(A6) 



Assuming CP conservation, this transformation is an orthogonal transformation. We shall denote the eigenstates of 
the neutralino mass matrix by XiJ X2; Xsj X4 with eigenvalues Al^o ^ ^ 3 4' labeled in order of increasing mass. Explicit 
expressions for these can be obtained, but these are not very illuminating. The neutralinos are mixtures of gauginos 
and Higgsinos 



(A7) 



One can obtain information on the neutralino masses by studying the expansion of the neutralino mass matrix (|A5[) 
in terms of Mz/ H for Mz <C This expansion is obtained most conveniently by using the basis (—17, —iZ'^, H'^, H^)i 
where 



1 



7 = 






ZO = 


1 










Hi = 


1 









Ht = 


1 









(,9'A3+.9A'), 
(.9A3-5'A'), 



(A8) 
(A9) 
(AlO) 
(All) 



are the photino, zino, and linear combinations of Higgsino states. In this basis, after a similarity transformation (see 
e.g. 3), the neutralino mass matrix can be written as 



/ Ah -AIz cos (/3 - f ) sw Mz sin (/3 - f ) \ 

M2 AIzCos{p-l)cw -AIzsin{p~l)cw 

-AIzcos{p~^)sw Mzcos{p-l)cw 

\ Mz sin - ^)sw -A/z sin (/3 - f ) civ -fi / 



(A12) 



The mass matrix (|A12[) can be diagonalized by using perturbation theory for values of AIz <C /i. For the case 
Ah < Ah, which is what one obtains in gravity mediated supersymmetry breaking (see below), the mass of the 
lightest neutralino can be written as, up to terms of 0{Mz / fJ.)^ , 



Ah 



sin 2 13 



1 



Similarly, for the second lightest neutralino X2 one obtains 



MhlMi 



M2 - Ml 



sin^ 2/3 



X2 



Mo 



Mlc 



sin 2/3 - -^j ( M^c^Ma 



Mfllfk sin^ 2/3 

Ml - A h 



(A13) 



(A14) 
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where = cos^ 9w and = sin^ 6iy. If instead we have M2 < Mi, a situation that arises in anomaly mediated 
supersymmetry breaking models, Eq. (|A14p would represent the mass of the lightest neutralino Xi^ aiid Eq. (|A13p 
would give the formula for the mass of the second lightest neutralino. The dependence of the lightest neutralino mass 
on the specific SUSY breaking scenario is due to the fact that the ordering of the gaugino mass parameters is model 
dependent. The leading term in the remaining two neutralino masses is the Higgs(ino) mixing parameter Thus if 
1^1 value is small compared to A/1,2, Higgsino can form a large or even dominant component of the lightest neutralino, 
as can be seen from the mass formulae for the remaining two ncutralinos: 



- + + ---^^) (,-Mi)(,-M2) + 8# (^1^) 



and 



= -A*- ^(1 -sin 2/3 4 5^— i - -4 cos^ 2/3. (A16) 
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